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ABSTRACT 

We consider a sample of 51 distant galaxy clusters at 0.15 <; z <; 0.9 (< z >~ 0.3), each cluster having 
at least 10 galaxies with available redshift in the literature. We select member galaxies, analyze the 
velocity dispersion profiles, and evaluate in a homogeneous way cluster velocity dispersions and virial 
masses. 

We apply the same procedures already recently applied on a sample of nearby clusters (z < 0.15, 
Girardi et al. 1998b) in order to properly analyze the possible dynamical evolution of galaxy clusters. 
We remark problems induced by the poor sampling and the small spatial extension of the sampled cluster 
region in the computation of velocity dispersion. 

We do not find any significant difference between nearby and distant clusters. In particular, we consider 
the galaxy spatial distribution, the shape of the velocity dispersion profile, and the relations between 
velocity dispersion and X-ray luminosity and temperature. Our results imply little dynamical evolution 
in the range of redshift spanned by our cluster sample, and suggest that the typical redshift of cluster 
formation is higher than that of the sample we analyze. 

Subject headings: galaxies: distances and redshifts - X-rays: galaxies - cosmology: observations. 



1 INTRODUCTION 

The knowledge of the properties of galaxy clusters and 
of their possible evolution plays an important role in the 
study of large scale structure formation constraining cos- 
mological models (e.g., Henry et al. 1992; Oukbir & Blan- 
chard 1992; Colafrancesco & Vittorio 1994; Eke, Cole, & 
Frenk 1996). The evolution of their statistical properties 
have been studied previously. In particular, there is no 
evidence of evolution in the bulk of population of X-ray 
selected clusters (out to z ~ 0.8; e.g., Burke et al. 1997; 
Jones et al. 1998; Rosati et al. 1998), with evidence for a 
negative evolution of the X-ray luminosity function hold- 
ing only for the brightest objects (e.g., Gioia et al. 1990; 
Vikhlinin et al. 1998; Rosati et al. 2000). There is ev- 
idence of a mild evolution of the cluster X-ray tempera- 
ture function (out to z ~ 0.8, Henry 1997; Donahue & 
Voit 1999), and of somewhat larger evolution of the in- 
ternal velocity-dispersion function (Carlberg et al. 1997b; 
Borgani et al. 1999). Other recent studies concern the 
relations between X-ray properties or between X-ray and 
optical properties finding no evidence of evolutions (out 
toz~ 0.4 - 0.5; e.g., Mushotzky & Scharf 1997; Borgani 
et al. 1999; Schindler 1999). Moreover, no evidence of 
evolution is found for other cluster properties, such as the 
iron abundance (out to z ~ 0.8 Mushotzky & Loewenstein 
1997) and the core radius of the distribution of hot intra- 
cluster medium (Vikhlinin et al. 1998). All such signs of 
evidence suggest a low value for the matter density param- 
eter n m (Carlberg et al. 1997b; Fan, Bahcall, & Cen 1997; 
Henry 1997; Borgani et al. 2000; see Mushotzky 2000 for 
a review). 

However, the validity of these studies relies on our actual 
understanding of the internal physics of both nearby and 
distant clusters. In particular, there is evidence that sev- 



eral clusters at moderate/distant redshift (z ~ 0.2 out to 
z ~ 1) are far from the state of dynamical equilibrium sug- 
gesting that present observations are reaching the epoch 
of cluster assembly. For instance, it is claimed that dis- 
tant clusters often show discrepancy in determination of 
mass estimates (e.g., Miralda-Escude & Babul 1995; Wu 
& Fang 1996; 1997), where the problems concern, in par- 
ticular, the cores of clusters (e.g., Allen 1998; Wu, Fang, 
Xue 1998), and are probably due to the lack of dynamical 
equilibrium or of spherical symmetry (e.g., Allen, Fabian, 
& Kneib 1996; Girardi et al. 1997b). Moreover, direct 
optical and X-ray observations show the strong elongation 
of some distant clusters (e.g., Gioia et al. 1999). 

In this framework, it is worth to analyze the internal 
dynamics of distant clusters comparing, in particular, the 
results with those obtained for nearby clusters. Here we 
focus our attention on the results as they come from the 
kinematical and spatial analysis of cluster member galax- 
ies. 

As for nearby clusters (at redshift z <; 0.15), avail- 
able results are based on very large samples, up to 100 
clusters, each with several galaxy redshifts available and 
treated in homogeneous way: the ENACS (ESO Nearby 
Abell Cluster Survey, Katgert et al. 1998) sample, and 
compilations collecting ENACS data and other clusters 
from the literature (den Hartog & Katgert 1996; Fadda 
et al. 1996, and the following updating by Girardi et al. 
1998b - hereafter F96 and G98, respectively). Significant 
substructures are found for 30-40% of clusters from both 
the distribution of member galaxies (e.g., Girardi et al. 
1997a; Biviano et al. 1997; Solanes, Salvador-Sol'e, and 
Gonzalez-Casado 1999) and X-ray analyses (Jones & For- 
man 1999), with a good one-to-one correspondence be- 
tween the optical and the X-ray images (Kolokotronis et 
al. 2000). However, with the exception of strongly 
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TABLE 1 
Cluster Sample 



Cluster Name 


Other Names 


N 


References 


(1) 


(2) 


(3) 


(4) 


A115 a 


ZwC10053. 4+2604 


28 


1 


A140" 


EDCC 520 


11 


2 


A222 




33 


3 


A223 




28 


3 


A370 




58 


4 


A520 


MS0451. 5+0250 


27 


3 


A521 




49 


5 


A665 


ZwC10826.1+6554 


41 


6 


A851 


CI 0939+47 


137 


4,7 


A1300 




95 


8 


A1689 




130 


9 


A2218 




53 


10 


A2390 


RXJ2153.6+1741 


325 


11 


A2744 


AC118 


76 


12,13 


A3639 




14 


14 


A3854 


C52 


41 


15 


A3888 


CL22315-3800 


98 


9 


A3889 




26 


16 


AS506 a 


CL0500-24 


29 


17 


AS910 


AC103 


88 


12,13 


AS1077 


AC 114 


103 


13,18 


CL0017-20 a 




26 


17 


CLJ0023+0423 


GHO0021+0406 


107 


19,20 


CL0024+16 


ZwC10024.0+1652 


134 


4,7 


CL0053-37 




22 


16 


CL0054-27 


J1888.16CL 


25 


4 


CL0303+17 


GHO0303+1706 


84 


4 


CL0412-65 




24 


4 


CL0949+44 


GHO0949+4408 


33 


7 


CL1447+26 




29 


4 


CL1601+42 


GHO1601+4259 


101 


4,7 


CLJ1604+4304 


GHO1602+4312 


95 


19,20 


F1637.23TL 




19 


21 


F1652.20CR 




20 


21 


J2175.15TR 




19 


21 


J2175.23C 




19 


21 


MS0015.9+1609 


CL0016+1609 HSTJ001831+16207 


111 


22 


MS0302.7+1658 


C10302+1658 


96 


23 


MS0302.5+1717 


C10302+1717 


43 


24 


MS0440.5+0204 a 




56 


25 


MS0451. 6-0305 




113 


22 


MS1008.1-1224 




109 


26 


MS1054.4-0321 a 




32 


27 


MS1224.7+2007 




54 


28 


MS1358.4+6245 


ZwC11358.1+6245 


281 


26 


MS1512.4+3647 




282 


28 


MS1621. 5+2640 




262 


23 


RXJ1716+67 




37 


29 


lE0657-56 a 


RASS1 069 


32 


30 


3C206 a 




15 


31 


3C295 




38 


4,7 



a Clusters for which spectral/morphological information is not available. 

(1) Zabludoff et al. 1990; (2) Collins et al. 1995; (3) Proust ct al. 2000; (4) Dressier ct al. 1999; (5) Maurogordato et al. 1999; (6) Oegerle ct al. 1991; 
(7) Dressier & Gunn 1992; (8) Lemonon et al. 1997; (9) Teague et al. 1990; (10) Lc Borgne et al. 1992; (11) Yee et al. 1996a; (12) Couch ct al. 1998; 
(13) Couch & Sharpies 1987; (14) Garilli et al. 1991; (15) Colless & Hewett 1987; (16) Cappi et al. 1998; (17) Infante et al. 1994; (18) Couch ct al. 1994; 
(19) Lubin ct al. 1998b; (20) Postman et al. 1998; (21) Bower ct al. 1997; (22) Ellingson et al. 1998; (23) Ellingson ct al. 1997; (24) Fabricant et al. 
1994; (25) Gioia ct al. 1998; (26) Yce ct al. 1998; (27) Tran ct al. 1999; (28) Abraham et al. 1998; (29) Gioia ct al. 1999; (30) Tucker et al. 1998; (31) 
Ellingson ct al. 1989. 

Cluster names: "A" for the catalog of Abell et al. (1989) and, in particular, "AS" for the supplementary southern clusters; "AC" and "C" for clusters 
used by Couch & Newell (1984), and by Colless & Hewett (1987), respectively, taken from the southern extension of the Abell catalog (Abell et al. 1989) 
in preparation at that time; "EDCC" for the Edinburgh-Durham Cluster Catalog (Lumsden et al. 1992); "F" and "J" for the catalog of Couch ct al. 
(1991); "GHO" for the catalog of Gunn et al. (1986); "HST" for the Hubble Space Telescope Medium Deep Survey cluster sample of Ostrander et al. 
(1998); "MS" for the Extended Medium Sensitivity Survey (Gioia et al. 1990); "RASS1" for the sample of bright clusters of galaxies in the southern 



substructured clusters (e.g., a 10% of bimodal clusters, see 
Girardi et al. 1997a; G98), most clusters seem not to be 
far from a global dynamical equilibrium. Finally, galaxy 
light is a good tracer of dark matter (e.g., Natarajan et 
al. 1998). The comparison between reliable estimates of 
velocity dispersions and X-ray temperatures of clusters 
suggests that the galaxy and hot gas components are not 
far from energy equipartition per unit mass; the possible 
discrepancies are likely to require extra-heating sources 
for poor clusters (e.g., White 1991; Bird, Mushotzky, & 
Mctzlcr 1995; G98). There is an overall agreement be- 
tween mass estimates inferred from the analysis of mem- 
ber galaxies and those from X-ray analysis of the hot gas 
(G98). Other detailed studies concern comparative analy- 
ses of different galaxy populations and their use as tracers 
of the cluster potential (e.g., early- and late-type ones, 
galaxies with or without emission lines, see, e.g., Biviano 
et al. 1997; Adami, Biviano, & Mazure 1998). 

As for more distant clusters (z > 0.2), the possible evo- 
lution of member galaxies is well studied (e.g., Butcher & 
Oemler 1978; Dressier et al. 1997; Abraham et al. 1998), 
but there are less definitive results on cluster internal dy- 
namics. Most results come from the analysis of the 16 
clusters at intermediate redshifts of CNOC (Canadian Net- 
work for Observational Cosmology, 0.18 < z < 0.55 ; Yee, 
Ellingson, & Carlberg 1996b) which represent a remark- 
ably homogeneous sample. In particular, as also found in 
nearby clusters, Lewis et al. (1999) claim for consistency 
between masses coming from optical and X-ray data, and 
Carlberg et al. (1997a) find that blue and red galaxies have 
different distributions in velocity and position. However, 
the difficulty of obtaining many redshifts in distant clusters 
has prevented from building larger samples. Rather, sev- 
eral works, concerning one or a small number of clusters, 
and using different techniques of analysis, can be found in 
the literature. 

The availability of a variety of techniques, already ap- 
plied to nearby clusters, suggests their application to dis- 
tant clusters. We thus ensure the homogeneity of our re- 
sults over a large range of cosmological distances. A ho- 
mogeneous analysis is in fact fundamental for the under- 
standing of the evolution of cluster properties. 

Here, we apply the techniques already used by G98 (cf. 
also F96) on a sample of 170 nearby clusters (at z < 0.15, 
data from ENACS and other literature) to analyze a col- 
lection of 51 distant clusters at 0.15 <; z <; 0.9. 

The paper is organized in the following manner. We 
shortly describe the data sample and our selection pro- 
cedure for cluster membership assignment in § 2 and § 3, 
respectively. We compute internal velocity dispersions and 
masses for clusters in § 4, with the exception of the three 
clusters with strong dynamical uncertainties which are dis- 
cussed in the Appendix A. We compare the "active" and 
the non "active" galaxies in § 5. We compare our results 
with those coming from X-ray and weak gravitational lens- 
ing analyses in § 6. We give a brief summary of our main 
results and draw our conclusions in § 7. 

Unless otherwise stated, we give errors at the 68% con- 



fidence level (hereafter c.l.) 

A Hubble constant of 100 h km s _1 Mpc^ 1 and a de- 
celeration parameter of qo = 0.5 are used throughout. 

2 THE DATA SAMPLE 

We consider 51 clusters at moderate/distant redshift 
z > 0.15 (median z = 0.33), each cluster having at least 10 
galaxies with available redshift and showing a significant 
peak in the redshift space. We remark that, due to the the 
difficulty of obtaining redshift data for distant clusters, we 
relax the requirements already applied to the sample of 
nearby clusters of G98 (z < 0.15, cf. also F96). In partic- 
ular, we consider also clusters with less than 30 available 
galaxy redshift, although we never take into account clus- 
ters which, after the procedure of the rejection of interlop- 
ers, are left with < 5 member galaxies. We relax also other 
requirements concerning the reliability of the estimate of 
velocity dispersion, i.e. the requirement of small errors on 
velocity dispersion (<; 150 km s" 1 ), and of flat integrated 
velocity dispersion profile in external cluster regions (cf. 
§4.1 for more details). 

Cluster data are collected from the literature. In or- 
der to achieve sufficiently homogeneous cluster data, the 
galaxy redshifts in each cluster are usually taken from 
only one reference source; different sources are used only 
when the data-sets are proved to be compatible. The 
data used for each cluster concern galaxy positions, red- 
shifts with the respective errors, and, when available, spec- 
tral/morphological information. 

Table 1 lists all the 51 clusters considered: in Col. (1) we 
list the cluster names; in Col. (2) we report other alterna- 
tive names found in the literature; in Col. (3) the number 
of galaxies with measured redshift in each cluster field; in 
Col. (4) the data references. 

3 CLUSTER MEMBER SELECTION 

In order to select member galaxies, we apply the same 
procedure as G98 (cf. also F96). 

Wc first use position and velocity information sequen- 
tially; then we use the two sets of data combined. In the 
first two steps we use the adaptive kernel technique by 
Pisani (1993,1996) as described in the Appendix A of Gi- 
rardi et al. (1996). The adaptive kernel technique is a 
nonparametric method for the evaluation of the density 
probability function underlying an observational discrete 
data set. For each detected peak, the method gives the 
corresponding significance and object density, as well as 
the associate objects. 

Firstly, we apply to each cluster field the two- 
dimensional adaptive kernel analysis to detect clusters 
which show an obvious bimodality in their projected 
galaxy distribution: i.e. formed by two significant (> 99% 
c.l.) clumps separated by a distance of ;> 0.5 h' 1 Mpc . 
These clumps are then analyzed separately. 

Afterwards, we apply the one-dimensional analysis to 
find the significant peaks in velocity distributions. The 
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TABLE 2 
Cluster Membership 



Cluster Name N p N g N m < z > Center 

R.A.(J2000) <5(J2000) 



(1) 


(2) 


(3) 


(4) 


(5) 








(6) 






A115S 


16 


13 


13 


.1958 


00: 


:55 


:59 


.90 


+26: 


:20: 


08.3 


A140 


7 


7 


7 


.1600 


01 


01 


:31 


55 


-23: 


:57. 


:46.2 


A222 


29 


27 


26 


.2138 


01: 


37 


:33 


83 


-12: 


:59 


:21.1 


A223 


18 


18 


14 


.2119 


01: 


:37: 


56 


.48 


-12: 


.48 


:33.9 


A370 


38 


37 


35 


.3744 


02 


:39 


:51 


.58 


-01: 


:34: 


:12.4 


A520 


19 


18 


18 


.2000 


04: 


:54: 


:14 


42 


+02: 


:57: 


:14.8 


A521 


40 


37 


35 


.2474 


04: 


:54: 


09 


.13 


-10: 


.11: 


:25.1 


A665 


31 


26 


25 


.1806 


08: 


30 


:46. 


85 


+65: 


53: 


:52.9 


A851 


67 


65 


55 


.4061 


09 


:42: 


:58 


.15 


+46: 


:59: 


:34.9 


A1300 


62 


59 


53 


.3078 


11 


:31: 


:57 


.69 


-19: 


:54: 


:35.2 


A1689a 


41 


38 


38 


.1837 


13: 


:11 


31 


62 


-01: 


:20: 


:58.0 


A1689b 


16 


16 


15 


.1746 


13: 


:11 


:28. 


58 


-01: 


:20: 


25.3 


A1689ab 


57 


50 


49 


.1821 


13: 


:11 


:30 


.24 


-01: 


:20: 


:54.2 


A2218 


45 


43 


43 


.1761 


16: 


:35 


:51 


96 


+66: 


:12: 


.19.8 


A2390 


243 


211 


200 


.2282 


21: 


:53 


:36 


80 


+ 17: 


41: 


:32.2 


A2744a 


36 


36 


34 


.3014 


00: 


:14: 


:21 


.26 


-30: 


.23 


:49.3 


A2744b 


27 


26 


25 


.3148 


00 


:14 


:20 


57 


-30: 


:24: 


:04.4 


A2744ab 


63 


57 


55 


.3078 


00: 


:14 


:21. 


16 


-30: 


:23: 


:52.4 


A3639 


8 


8 


7 


.1480 


19 


:28 


:18 


.41 


-50: 


:54: 


:28.6 


A3854a 


21 


18 


18 


.1520 


22: 


17 


:53 


03 


-35: 


:46: 


:42.7 


A3854b 


13 


13 


9 


.1459 


22: 


■17 


:33 


97 


-35: 


:44 


:47.1 


A3854ab 


34 


34 


30 


.1506 


22: 


:17: 


.46 


.50 


-35: 


:45 


:12.3 


A3888 


81 


55 


50 


.1508 


22: 


:34: 


:26. 


.90 


-37: 


:43: 


:51.0 


A3889a 


7 


7 


7 


.2559 


22: 


:34: 


:54. 


.47 


-30: 


33 


:50.5 


A3889b 


9 


9 


9 


.2495 


22 


:34 


:49. 


.35 


-30: 


:32 


:13.1 


AS506 


23 


21 


21 


.3201 


05: 


01 


:11 


.85 


-24 


:25: 


:01.5 


AS910 


56 


54 


53 


.3076 


20 


:57 


02 


89 


-64 


:40: 


:04.7 


AS1077 


85 


70 


63 


.3148 


22: 


:58 


:47. 


.14 


-34: 


.47 


:59.8 


CL0017-20 


20 


20 


20 


.2717 


00: 


:19 


37 


.82 


-20: 


:26: 


39.1 


CLJ0023+0423a 


14 


14 


5 


.8453 


00: 


:23: 


52 


.69 


+04: 


19: 


:38.3 


CLJ0023+0423b 


7 


7 


5 


.8273 


00 


:23: 


53 


.82 


+04: 


•23: 


:16.2 


CL0024+16 


102 


95 


73 


.3937 


00: 


:26 


:34. 


79 


+ 17: 


10 


.04.8 


CL0053-37 


20 


20 


20 


.1652 


00: 


:55: 


:59 


.44 


-37: 


32 


:36.1 


CL0054-27 


10 


9 


7 


.5604 


00: 


:56: 


56 


.04 


-27: 


.40 


:31.9 


CL0303+17 


44 


43 


29 


.4195 


03: 


:06 


11 


13 


+ 17: 


18: 


:09.0 


CL0412-65 


7 


7 


6 


.5086 


04: 


:12: 


:53 


.39 


-65: 


:51 


:13.2 


CL0949+44a 


15 


15 


14 


.3781 


09: 


:52: 


:57 


.50 


+43: 


:55: 


:37.0 


CL0949+44b 


9 


9 


8 


.3493 


09: 


:53 


:01 


25 


+43: 


:55: 


22.6 


CL1447+26 


16 


11 


5 


.3763 


14: 


:49 


.28 


.78 


+26: 


06: 


58.3 


CL1601+42 


57 


53 


46 


.5403 


16: 


:03: 


:10 


.46 


+42: 


:45: 


:37.2 


CLJ1604+4304 


19 


14 


8 


.9018 


16: 


:04 


:25. 


.09 


+43: 


04: 


:11.0 


F1637.23TL 


8 


8 


6 


.47903 


23: 


:59: 


:20 


.68 


-32: 


:17: 


:45.3 


F1652.20CR 


8 


8 


6 


.4102 


04: 


:47: 


:57 


.66 


-20: 


:37: 


:29.6 


J2175.15TR 


8 


8 


8 


.3948 


03: 


:34: 


:20 


71 


-38: 


:53: 


:54.5 


J2175.23C 


10 


6 


5 


.4058 


03: 


:32: 


59 


.30 


-39: 


.06: 


:49.7 


MS0015.9+1609 


63 


50 


42 


.5490 


00: 


:18 


33 


.49 


+ 16: 


26: 


:02.5 


MS0302. 7+1658 


37 


33 


30 


.4248 


03: 


:05 


31 


.63 


+ 17: 


10: 


.12.0 


MS0302. 5+1717 


28 


26 


24 


.4241 


03: 


:05: 


:17 


.92 


+ 17: 


'28: 


:34.4 


MS0440.5+0204 


32 


32 


32 


.1969 


04: 


:43: 


09 


.46 


+02: 


10: 


29.5 


MS0451. 6-0305 


67 


46 


40 


.5403 


04: 


:54: 


:11 


24 


-03: 


:00 


:45.4 


MS1008. 1-1224 


74 


71 


65 


.3070 


10: 


:10: 


:31 


.76 


-12: 


.40: 


:05.4 


MS1054.4-0321 


32 


32 


32 


.8318 


10 


:56: 


:57 


.31 


-03: 


37 


:44.2 


MS1224.7+2007 


23 


23 


23 


.3253 


12: 


:27 


:18. 


.81 


+ 19: 


■50: 


:26.7 


MS1358.4+6245 


185 


141 a 


133 


.3278 


13 


:59: 


:50 


.92 


+62: 


:30: 


:49.8 


MS1512.4+3647 


70 


46 


35 


.3711 


15: 


:14: 


:16 


.45 


+36: 


:34: 


:57.7 


MS1621.5+2640 


119 


106 


88 


.4271 


16: 


23 


:34. 


.52 


+26: 


:34: 


:17.1 


RXJ1716+67 


37 


32 


19 


.8073 


17: 


:16: 


:48 


.92 


+67: 


08: 


:21.1 


1E0657-56 


14 


12 


12 


.2966 


06: 


:58: 


:37 


.83 


-55: 


:56 


:56.0 


3C206 


7 


A 


7 


.1980 


08: 


39 


50 


.10 


-12: 


11 


32.2 


3C295 


21 


15 


.4591 


14: 


:11 


:21 


.54 


+52: 


:11: 


54.6 



main cluster body is naturally identified as the highest sig- 
nificant peak. All galaxies not belonging to this peak are 
rejected as non cluster members. F96 and G98 required 
that peaks are significant at the 99% c.l. and, for clusters 
with secondary peaks, they assumed that the peaks are 
separable when their overlapping is < 20% and their ve- 
locity separation is Aw > 1000 km (here, we consider 
1000 km s _1 in the appropriate cluster rest-frame). In 
dealing with distant clusters, we apply the peak analysis 
to very poor samples, so obtaining small peak probability: 
in a few fields we identify clusters with the highest peak 
having significance < 99% (but always > 95%). 

The combination of position and velocity information, 
represented by plots of velocity vs. clustercentric distance, 
reveals the presence of surviving interlopers (e.g., Kent & 
Gunn 1982; Regos & Geller 1989). To identify these inter- 
lopers in the above-detected systems we apply the proce- 
dure of the "shifting gapper" by F96. We apply the fixed 
gap method to a bin shifting along the distance from the 
cluster center. According to F96 prescriptions, we use a 
gap of > 1000 km s _1 (in the cluster rest-frame) and a 
bin of 0.4 hr 1 Mpc, or large enough to include 15 galaxies. 
As for very poor distant clusters (with less than 15 mem- 
bers), we reject galaxies that are too far in velocity from 
the main body of galaxies of the whole cluster (considering 
a somewhat larger gap, cf. § 3.1 for details). 

When early- and late-galaxy type populations showed 
different mean and variance in the velocity distribution, 
Girardi et al. (1996) retained only the early population 
as better tracer of the cluster potential. Similarly, when 
only spectral information was available, F96 applied the 
same procedure by rejecting emission line galaxies (here- 
after ELGs, while NELGs indicate galaxies without emis- 
sion lines). Indeed there is evidence that ELGs lead to too 
high estimates of internal velocity dispersion (e.g., Koranyi 
& Geller 2000) and that they could be not in dynami- 
cal equlibrium within the cluster potential (Biviano et al. 
1997). For distant clusters, galaxy morphologies are gen- 
erally not available, but spectral type and/or color give 
information about the presence of some nuclear galaxy 
activity or strong star formation. For the 43 out of 51 
cluster samples with available information, we reject the 
likely "active" galaxies (hereafter AGs), i.e. galaxies hav- 
ing a strong star formation activity and/or signs of nuclear 
activity (see our classification in § 3.2). 

Table 2 lists the results of the member selection proce- 
dure (cf. § 3.1 for details on some specific clusters). In 
Col. (1) we list the system name, i.e. the name of the 
parent cluster with possible indication of the peak (e.g. 
A1689a and Al689b); in Col. (2) the number of galaxies 
found by the adaptive kernel method in each peak, N p ; 
in Col. (3) the number of galaxies left after the "shifting 
gapper", N g ; in Col. (4) the number of member galaxies 
after the rejection of the AGs, N m , and used to compute 
the mean redshift determined via the biweight estimator 
(Beers, Flynn, & Gebhardt 1990), and the cluster center 
as determined via the two-dimensional adaptive kernel (in 
Cols. (5) and (6), respectively). 



3. 1 Results for Specific Clusters 

According to our analysis of galaxy distribution, we find 
indication for bimodality only in the case of A115 (cf. also 
Forman et al. 1981; Beers, Geller, & Huchra 1983). More- 
over, we consider only the Southern peak, A115S, since 
the Northern peak has too few galaxies to survive to the 
whole procedure of member selection. 

According to our analysis of cluster velocity distribu- 
tions, out of the 51 cluster fields here analyzed, we find 
51 well-separated peaks (45 from 45 one-peaked fields and 
six from three two-peaked fields) and three fields with two 
strongly superimposed peaks. 

The three fields showing two separable peaks are 
A3889, CL0949+44 (cf. also Dressier & Gunn 1992), 
and CLJ0023+0423. In particular, the field of the 
CLJ0023+0423 cluster shows a complex structure in the 
velocity distribution, containing a system of four peaks 
strongly superimposed; however, when reanalyzing only 
galaxies belonging to this system, we find two peaks. 
These peaks corresponds to those found by Lubin, Post- 
man, & Oke (1998a, CLJ0023+0423 "A" and "B" instead 
of our "b" and "a"). 

The cluster fields for which the peak separation is not 
secured at a high c.l. arc: A1689 (cf. also Girardi et al. 

1997) , A2744, and A3854. The strongly overlapped peaks 
could indicate the presence of substructures in a single 
system and, in this case, the dynamics of these clusters 
is strongly uncertain; therefore we consider both the case 
with the two peaks disjoined or together (e.g., A1689a, 
A1689b, and Al689ab), cf. the Appendix A for other de- 
tails. 

As for very poor distant clusters (with less than 15 
galaxies), the procedure of the "shifting gapper", which 
works with a gap of 1000 km s _1 in a shifting bin con- 
sidering at least 15 galaxies, cannot be applied. In these 
cases we reject galaxies that are too far in velocity from 
the main body of galaxies of the whole cluster rather than 
in a shifting bin. Moreover, we adopt a slightly larger 
gap since the suitable size of the gap increases with the 
available statistics (cf. the density-gap by Adami et al. 

1998) . We reject one galaxy in CL0054-27, one galaxy 
in J2175.23C, and two galaxies in 1E0657-56, where the 
gap is ;> 2000 km s^ 1 . The situation is less obvious for 
other two galaxies in CL0054-27 and other three galax- 
ies in J2175.23C, where the respective gaps are ~ 1150 
km s" 1 , - 1300 km s" 1 , respectively. For CL0054-27, the 
two uncertain members are close to the center and we de- 
cide to retain them. Instead, for J2175.23C, we decide to 
reject the three uncertain members, which are connected 
to the main body of galaxies only thanks to the presence of 
an AG. Moreover, two of the uncertain members are AGs, 
and the other uncertain member would result the most 
distant galaxy from the cluster center. Monte Carlo sim- 
ulations performed in § 4.4 show that, in the case of poor 
statistics, a fixed gap of ~ 1250 km s^ 1 allows us to well 
recover, on average, the estimate of velocity dispersion. 

As for the combined analysis of position and velocity in- 
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formation, the plot of rest-frame velocity versus projected 
clustercentric distance of MS1358. 4+6245 shows the exis- 
tence of a close system corresponding to a southern group 
(cf. Carlberg et al. 1996). We exclude this group reject- 
ing all galaxies outside 1.2 h^ 1 Mpc (cf. also Borgani et 
al. 1999). 

3.2 Classification of "Active Galaxies" 

As for our classification of AGs, in several cluster sam- 
ples only main galaxy spectral features are reported. In 
these cases we classify as "active" galaxies to be rejected 
those where the presence of emission lines is reported. For 
other clusters, where more detailed information is given, 
we always reject galaxies with very strong emission lines or 
classified as starbust or AGN: in the following we describe 
the classification adopted for these specific studies. 

As for the data by Postman, Lubin, & Oke (1998), we 
reject galaxies with the presence of a [Oil] line with an 
equivalent width of EW[OII] ;> 15 A , which corre- 
sponds, according to the authors, to an active, star forming 
galaxy. 

As for the data by Dressier et al. (1999), we con- 
sider as AGs those galaxies classified from their spectra 
as "e(a)" (with strong Balmer absorption plus [Oil] emis- 
sion), "e(n)" (having AGN spectra), or "e(b)" (with very 
strong [Oil] emission, possibly starbust galaxies). As for 
data by Dressier & Gunn (1992), "active" galaxies are 
those classified with "e" (with emission lines, usually [Oil] 
or [OIII]) or "n" (with very strong emission, likely due to 
AGN). 

As for the CNOC clusters (Yee et al. 1996a; Ellingson 
et al. 1997; Abraham et al. 1998; Ellingson et al. 1998; 
Yee et al. 1998), we consider as AGs those galaxies classi- 
fied from their spectra with "5" (with emission lines, likely 
irregular galaxies) or "6" (likely AGN/QSO). 

Finally, we also consider as AGs those galaxies labeled 
"starbust" by Couch et al. (1998, and previous works by 
Couch & Sharpies 1987 and Couch et al. 1994), which are 
characterized by blue colors and emission-filled US lines. 

4 ANALYSIS OF 
CLUSTER INTERNAL DYNAMICS 

The analysis of the three clusters which show strongly 
superimposed peaks (A1689, A2744, and A3854) is post- 
poned to the Appendix A. In the following sections we 
analyze the 45 cluster fields which show only one peak 
in their velocity distribution and the three cluster fields 
which show two separable peaks for a total of 51 well 
defined systems (cf. § 3.1). This sample can be compared 
to that of 160 well-separated peaks for nearby clusters of 
G98. We remark that, on average, the distant clusters are 
less well sampled as for both the number of cluster mem- 
bers (median N m = 21 vs. 39), and the spatial extension 
Rmax, which is the clustercentric distance of the most dis- 
tant galaxy from the cluster center (median R m ax — 0.64 
vs. 1.45 h" 1 Mpc). Throughout this analysis we apply 
homogeneous procedures already used by G98 (cf. also 



F96). 

4-1 Velocity Dispersions 

We estimate the "robust" velocity dispersion line-of- 
sight, cr v , by using the biweight and the gapper estimators 
when the galaxy number is larger and smaller than 15, re- 
spectively (cf. ROSTAT routines - see Beers et al. 1990), 
and applying the relativistic correction and the usual cor- 
rection for velocity errors (Danese, De Zotti, & di Tullio 
1980). In particular, for a few cases where the velocity 
error is not available, we assume a typical velocity error of 
300 km s _1 . When the correction for velocity errors leads 
to a negative value of a v , we list a v = 0. 

Following F96 (cf. also Girardi et al. 1996) we an- 
alyze the "integral" velocity dispersion profile (hereafter 
VDP), where the dispersion at a given (projected) radius 
is evaluated by using all the galaxies within that radius, 
i.e. cr v (< R). The VDPs allow to check the robustness of 
cr v estimate. In particular, although the presence of ve- 
locity anisotropy in galaxy orbits can strongly influence 
the value of a v computed for the central cluster region, 
it does not affect the value of the a v computed for the 
whole cluster (e.g., Merritt 1988). The VDPs of nearby 
clusters show strongly increasing or decreasing behaviors 
in the central cluster regions, but they are flattening out 
in the external regions (beyond ~1 hr 1 Mpc, cf. also den 
Hartog & Katgert 1996) suggesting that in such regions 
they are no longer affected by velocity anisotropics. Thus, 
while the ovvalues computed for the central cluster re- 
gion could be a very poor estimate of the depth of clus- 
ter potential wells, one can reasonably adopt the o~ v value 
computed by taking all the galaxies within the radius at 
which the VDP becomes roughly constant. As for the dis- 
tant clusters we analyze, when the data are good enough, 
the VDPs show a behavior similar to that of nearby clus- 
ters (cf. Figure 1). Unfortunately, distant clusters suffer 
for the poor sampling, and also for the small spatial ex- 
tension of the sampled cluster region. Indeed, the strongly 
decreasing VDP in the external sampled regions of some 
clusters (maybe the striking cases are AS506, CL0017-20, 
CL0054-27, 3C295) suggests that the correct estimates of 
velocity dispersions could be smaller than those, o~ v , we 
can estimate with present data; therefore, in these cases, 
o~ v should be better interpreted as an upper limit (see also 
some cases in F96; these cases were then not considered in 
G98). In other cases, when the member galaxies are too 
few, the analysis of VDPs does not allow any conclusion. 

In Table 3 we report the value of a v computed consid- 
ering all member galaxies. However, we make a note for 
the clusters which do not share the requirements for the 
nearby clusters of F96 and G98, i.e.: with an original num- 
ber of galaxies in the field smaller than 30; with a peak in 
the velocity distribution less significant than 99%; with an 
error on a v larger than 150 km s , with a VDP which is 
poorly defined or without a flat behavior in the external 
cluster regions. 

After fixing the cosmological background, the theory of 
a spherical model for nonlinear collapse allows to recover 
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the value of the radius of virialization, R V i r , within which 
the cluster can be considered not far from a status of dy- 
namical equilibrium. The relation between the density of 
a collapsed (virialized) region and the cosmological den- 
sity is p vir — \&n 2 p cr — 187r 2 x 3H 2 /8itG (for a Q m = 1 
universe). As a first approximation, the mass contained 
within R v i r , M V i r = (47r/3) ■ R 3 ir p V i r , is given by the virial 
estimate (3n/2) ■ (a 2 R vir /s / 'G) , where /s depends on the 
details of galaxy spatial distribution (e.g., G98). There- 
fore R 2 mr = alh/(6TrH 2 ), where H = 100 • h(l + zfl 2 . 
For nearby clusters G98 give a first roughly estimate of 
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Fig. 1. — Integrated line-of-sight velocity dispersion profiles 
a v {< R), where the dispersion at a given (projected) radius 
from the cluster center is estimated by considering all galax- 
ies within that radius. The bootstrap error bands at the 68% 
c.l. are shown. The horizontal lines represent X-ray tempera- 
ture with the respective errors (cf. Table 4) transformed in a v 
imposing f3 spec — 1 (where f3 spec = a% / (kT / pm p ) , with p, the 
mean molecular weight and m p the proton mass). The vertical 
faint lines indicate the virialized region within R v i r . 

Rvir ~ 0.002 • a v (km~ 1 s h^ 1 Mpc). A following re- 
estimate of Girardi et al. (1998a) suggests rather a scaling 
factor of 0.0017. Since we find that distant clusters have 
a galaxy distribution similar to that of nearby ones (see 
in the following), we adopt here the same scaling relation 
with a v : 

R mr - 0.0017 • a v /(l + zf 12 (km^s IT 1 Mpc) (1) 

introducing only the scaling with redshift (cf. also Carl- 
berg et al. 1997c for a similar relation). 

4-2 Galaxy Distribution 

As for the study of the spatial distribution of galaxies 
within distant clusters, following G98 (cf. also Girardi et 
al. 1995; Adami et al. 1998) we fit the galaxy surface 
density of each cluster to a King-like distribution (com- 
parable to the /3-profile in fitting the distribution of hot 
diffuse intracluster gas): 



[1 + (R/R c ) 2 ] a ' 



(2) 



where R c is the core radius and a is the parameter which 
describes the galaxy distribution in external regions (a = 1 
corresponds to the classical King distribution). This sur- 
face density profile corresponds to a galaxy volume-density 
p = p /[l + (r/R c ) 2 } 3 ^"' /2 , with f3 fu , 9 ai = {2a + l)/3, 
i.e. p(r) cx r~ 3l3fit - sal for r >> Rc- We perform the fit 
through the Maximum Likelihood technique, allowing Rc 
and a to vary from 0.01 to 1 h^ 1 Mpc and from 0.5 to 1.5, 
respectively. 

To avoid possible effects due to the non circular sam- 
pling of clusters, by visual inspection of the original sam- 
pled region of each cluster we extract the largest circular 
region, with center as in Table 2, there inscribed. We per- 
form the fit within this circular cluster region whose radius 
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we define as R m ax,c- We consider only the 30 clusters with 
at least ten member galaxies within R ma x,c and, in partic- 
ular, a subsamplc of 13 clusters with R max ,c/ Rvir > 0.5. 

The median value of a, with the respective errors at 
the 90% c.L, is = 0.63±g;g|, and 0.67 is found when we 
consider only the 13 clusters with a large sampled radius. 



This value agrees with a = 0.70_ 



0.03 



found for nearby 



clusters, and corresponds to a Pfu, ga l ~ 0.8, i.e. to a vol- 



ume galaxy-density p oc 



After fixing a — 0.7, we 



again fit the galaxy distribution of each cluster, obtaining 



a median value of R r = 0.045 



+0.005 
-0.015 



h- 1 Mpc (and 0.05 



h^ 1 Mpc for the well-sampled 13 clusters). Thus, in our 
cluster sample, the typical value of R c (and R V i r /R c ~ 20) 
is again in agreement with that found in nearby clusters 
where R c = 0.05 ± 0.01 h^ 1 Mpc. Hereafter, we assume 
the above King-modified distribution, with the same pa- 
rameters of nearby clusters, i.e. a — 0.7 and R c — 0.05 
h^ 1 Mpc, for all clusters of our sample. 

4-3 Virial Masses and Velocity Anisotropics 

Assuming that clusters are spherical, non rotating sys- 
tems, and that the internal mass distribution follows 
galaxy distribution, cluster masses can be computed 
throughout the virial theorem (e.g., Limber & Mathiews 
1960; The & White 1986) as: 



M = M v -C = — ■ 
where the projected virial radius, 



3-7T a^Rpv 



G 



-c, 



R PV = N(N-l)/(E i>j RT j 1 ), 



(3) 



(4) 



describes the galaxy distribution and is computed from 
projected mutual galaxy distances, Rij\ C is the surface 
term correction to the standard virial mass My and it is 
due to the fact that the system is not entirely enclosed in 
the observational sample (cf. also Carlberg et al. 1996; 
G98). 

Following G98 we want to estimate cluster masses con- 
tained within the radius of virialization, R V i r . In fact, 
clusters cannot be assumed in dynamical equilibrium out- 
side R V i r and considering small cluster regions leads to 
unreliable measures of the potential (o~ v could be strongly 
affected by velocity anisotropies) and of the surface term 
correction (Koranyi & Gcllcr 2000). 

Only few distant clusters are sampled out to R V i r . As for 
a v , the above analysis of the VDP give indications about 
its reliability, i.e. VDPs which are flat in the external 
cluster regions will give reliable estimates of a v even when 
clusters are not sampled out to R V i r . As for Rpy, which 
describes the galaxy spatial distribution, it can be recov- 
ered in an alternative theoretical way from the knowledge 
of the parameters of the King-like distribution (Girardi et 
al. 1996; see also G98 for a simple analytical approxima- 
tion in the case of a = 0.7 and R c /R V i r = 0.05). This pro- 
cedure allows to compute Rpy at each cluster radius and, 
in particular, we compute Rpv at R v i r , which is needed in 
the computation of the mass within R V i r . By using well- 
sampled nearby clusters G98 verified the reliability of this 



alternative estimate and evaluated the typical error intro- 
duced by the use of the average King-like parameters of 
the sample (~ 0.2 h^ 1 Mpc , corresponding to about 25% 
of R PV ). 

The computation of the C correction at the boundary 
radius, b, 



C = M V - 47Tb 3 



P(b) 



J Q 4irr 2 pdr 



[a r {b)/a{< b)f 



(5) 



requires the knowledge of the velocity anisotropy of galaxy 
orbits. In fact, a r (b) is the radial component of the veloc- 
ity dispersion a(b), while er(< b) refers to the integrated 
velocity dispersion within 6; here b — R V i r . Having as- 
sumed that in clusters the mass distribution follows the 
galaxy distribution, one can use the Jeans equation to es- 
timate velocity anisotropies from the data, i.e. from the 
(differential) profile of the line-of-sight velocity dispersion, 
a v (R). Unfortunately, this profile requires a large number 
of galaxies and we can compute it only by combining to- 
gether the data of many clusters, without preserving clus- 
ter individuality. 
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Fig. 2. — The (normalized) line-of-sight velocity dispersion, 
<j v (R), as a function of the (normalized) projected distance 
from the cluster center. The points represent data combined 
from all clusters and binned in equispatial intervals. We give 
the robust estimates of velocity dispersion and the respec- 
tive bootstrap errors. We give the results for distant clusters 
(open circles) and for nearby clusters taken from Girardi et 
al. (1998b, filled circles). The solid and dotted line represent 
the models for isotropic and moderate radial orbits of galaxies, 
respectively (see text). 

For both nearby and distant clusters, Figure 2 shows 
the observational a v (R) computed by combining together 
the galaxies of all clusters, i.e. by normalizing distances to 
R v i r and velocities, relative to the mean cluster velocity, 
to the observed global velocity dispersion a v . For nearby 
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TABLE 3 

Dynamical Properties 
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a Clusters having in their field less than 30 galaxies with available redshift (cf. Table 1). 
b Clusters with a peak in the velocity distribution less significant that 99%. 
c Clusters with an error on a v of > 150 km s -1 . 
d Clusters with a VDP which is poorly defined. 

e Clusters with a VDP which is without a flat behavior in the external cluster regions: the strong decreasing suggests that the estimates of a. 
My, and M are better interpreted as upper limits (see text). 9 



clusters the observational profile is well described by a the- 
oretical profile obtained by the Jeans equation, assuming 
that velocities are isotropic, i.e. that the tangential and 
radial components of velocity dispersion are equal (i.e., the 
velocity anisotropy parameter A = 1 — 0g(r) / '<J 2 r {r) = 0). 
For distant clusters this model is less satisfactory, but can- 
not be rejected being acceptable at the ~ 15% c.l. (accord- 
ing to the x 2 probability). 

In order to give C-corrections more appropriate to each 
individual cluster G98 used a profile indicator, I pi which 
is the ratio between a v (< 0.2 x R V i r ), the line of sight 
velocity dispersion computed by considering the galaxies 
within the central cluster region of radius R = 0.2 x R V i r , 
and the global a v . According to the values of this param- 
eter, they divided clusters in three classes containing the 
same number of clusters: "A" clusters with a decreasing 
profile (I p > 1.16), "C" clusters with a increasing profile 
(I p < 0.97), and an intermediate class "B" of clusters with 
very flat profiles (0.97 < I p < 1.16). Each of the three 
types of profiles can be explained by models with a differ- 
ent kind of anisotropy, i.e. radial, isotropic, and circular 
orbits in the case of A, B, and C clusters, respectively, re- 
quiring different values of [a r (R V i r )/a(< R V i r )] 2 , cf. G98. 
We can define 14, 11, and 8 clusters of class A, B, and 
C, respectively, and for each class we use the respective 
{<j r (R V i r )/<T(< R V ir)] 2 given by G98 to determine the C- 
corrections. The median values of the relative corrections 
C/M v are 45%, 20%, and 14% for A, B, and C clusters, 
respectively. For 18 clusters we cannot define the kind 
of profile and we assume the intermediate one. The me- 
dian correction on the whole sample is then C/My ~ 21%, 
very similar to that found by G98 for nearby clusters and 
to that suggested by Carlberg et al. (1997c) for CNOC 
clusters. 

In Table 3 we list the results of the cluster dynamical 
analysis: the number of cluster members as taken from Ta- 
ble 2, N m (Col. 2); the clustercentric distance of the most 
distant galaxy from the cluster center, R max (Col. 3); the 
global line-of-sight velocity dispersion a v with the respec- 
tive bootstrap errors (Col. 4); the radius which defines the 
region of virialization, R vir (Col. 5); the projected virial 
radius, Rpv, computed at R V i r (Col. 6); the cluster type 
according to their velocity dispersion profile, T (Col. 7); 
the estimate of cluster mass contained within R vir as de- 
termined from the standard virial theorem, My , and after 
the pressure surface term correction, M (Cols. 8 and 9, 
respectively). The errors on My take into account of the 
errors on a v and the above quoted error of 25% on Rpy. 
The percent errors on Mqv are the same as for My, i.e. 
we neglect uncertainties on C-correction. 

There are some possible suggestions that galaxy orbits 
in distant clusters have a somewhat more radial veloc- 
ity anisotropy than those in nearby clusters. Figure 2 
shows that the velocity dispersion profile for distant clus- 
ters seems less flat than that of nearby clusters. Indeed, 
according to the available data for distant clusters, also 
models with moderate radial orbits, i.e. A = r 2 /(r 2 + r 2 ) 
with r a = 0.25 x R v i r , can be acceptable (at the ~ 4% 



c.l.). Moreover, the distant clusters with strongly decreas- 
ing profiles (type A) are many more than clusters with 
increasing profiles (type C), while, for definition, the num- 
ber of A and C types are equal among nearby clusters. 
However, we verify that there is no significant difference 
between the combined profiles of distant and nearby clus- 
ters, and that distant and nearby clusters are not differ- 
ent according to the median value of the profile parame- 
ter I p . Therefore, with present data , we conclude that 
the possible evidence of a larger amount of radial velocity 
anisotropies in galaxy orbits of more distant clusters is not 
significant. 

4-4 Robustness of the Results 

Here we address the effects of the poor sampling on our 
results, in particular on our estimates of velocity disper- 
sions. 

In order to check the effects of the small number of red- 
shifts we perform Monte Carlo simulations by randomly 
undcrsampling the 20 cluster fields having more than 50 
galaxies and only one peak in the velocity distribution. 
Since our random undersampling does not consider other 
parameters like the proximity to the cluster center and the 
galaxy color, the following results should be considered as 
conservative. 

For each of the 20 cluster fields we perform 500 random 
simulations extracting ten galaxies each time (the lowest 
limit in our sample). Then we apply the whole procedure 
of member selection, i.e. the detection of a significant peak 
in the velocity distribution via the adaptive kernel method, 
the application of the fixed gap to the remaining galaxies, 
and finally the rejection of "active" galaxies. We use a 
fixed gap of 1250 km s _1 which seemed appropriate when 
a small number of galaxies is considered (cf. § 3.1). Only 
a fraction (~ 20%) of the simulated clusters survive to the 
procedure of the rejection of interlopers. These clusters 
contain typically 5-6 members to compute the velocity 
dispersion, a v s, and the associate statistical error, Aa v ^. 
For each of the 20 clusters we compute the median value, 
< o- v _i > ran , and the s.d., o~ ran {o~v,i), of the velocity dis- 
persions of the corresponding simulated clusters, as well 
as the median statistical error < Aa Vt i > rcm - 

We verify the robustness of our estimates of velocity 
dispersions by computing the median value and 90% c.l. 
errors of a v j < o~ v > ran within the sample of 20 clusters: 
we obtain < (a v / < a v > ran ) >20= l-Ollo.09 ( nxe d gaps 
of 1000 or 1500 km s' 1 give values of < (a v / < a v > ran 
) >2o= 1-10 or 0.97, again consistent with unity). 

These simulations also allow us to estimate the global er- 
ror on the estimate of the velocity dispersion, i.e. the error 
associated to the procedure of member selection in addi- 
tion to the statistical error connected to the selected mem- 
bers. On the whole sample of 20 clusters, the median value 
of the s.d. of velocity dispersions of simulated clusters is 
very high, < a ran (a v _i) >2o~ 850km s _1 , much larger 
than the corresponding statistical error, < (< Aa v s > ran 
) >20^ 350 km s _1 . Therefore, in the case of a very small 
number of available redshifts, the error associated to the 
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member selection procedure can be more important than 
the statistical one. However, this kind of error rapidly 
decreases as the amount of available data increases. For 
instance, we obtain a global error of ~ 500fcm s _1 vs. a 
statistical error of 300fcm s" 1 for simulated clusters of 15 
galaxies (on average 6-7 members) and the two error esti- 
mates become comparable for those clusters containing at 
least ten members. 

As for the poor spatial extension, as suggested by the 
VDPs of Figure 1 , the effect of individual clusters may be 
large. To be more quantitative, we consider the 11 clus- 
ters sampled out to R v i r and the corresponding velocity 
dispersions within R V i r /2: the estimate of velocity disper- 
sion varies by ~ 25% for three clusters. When considering 
the velocity dispersions within R V i r /A the variation is more 
than 25% for six clusters and reaches 65% for one cluster. 
Unfortunately, since the shape of VDPs range different 
behaviors, the effect cannot be predicted for individual 
clusters, although strongly decreasing/increasing profiles 
in the external sampled regions suggest that more cor- 
rect estimates of velocity dispersions would need data over 
larger field of view (cf., e.g. AS506, CL0017-20, CL0054- 
27, and 3C295). 

Another effect concerns the sampling within non circu- 
lar apertures. Since in our sample the elongation of the 
sampled region is not extreme (R m ax,c/ Rmax ~ 0.6), we 
expect that this effect is smaller than the previous one. 
One can quantify the effect by increasing the weight for ex- 
ternal cluster regions in the standard estimate of velocity 
dispersion a 2 me = [a 2 m (N t -l) + a 2 ve {N e -l)]/(N t + N e -l), 
where a V i and a ve are the velocity dispersions as computed 
on the internal and external cluster regions containing Ni 
and N e galaxies, respectively. Possible undersampling in 
the external regions can be corrected for by artificially in- 
creasing N e . We compute a V i and a ve inside and outside 
Rmax,c for the 35 clusters with at least 10 galaxies: a rea- 
sonable increase of the weight (by a factor four) for the 
external region leads to variations for individual clusters 
of the order of <; 7%. 

Finally, we check the effect of changing the classification 
of "active" galaxies. We consider the 22 clusters where au- 
thors classified also galaxies with ongoing moderate star 
formation, possibly spiral-like galaxies: for the data by 
Dressier et al. (1999) we consider "e(c)" galaxies (with 
moderate absorption plus emission, spiral-like); for the 
CNOC clusters we consider galaxies classified with "4" 
(Sbc), and for the clusters by Couch and collaborators we 
consider "Sp" galaxies (spiral-like, with spectra and color 
properties of normal nearby spiral galaxies). For these 22 
clusters we compare the velocity dispersions as computed 
in § 4.1, <t„, with those computed rejecting also spiral-like 
galaxies as defined above, a v -s- This different definition 
of "active" galaxies does not affect the estimate of velocity 
dispersions. We find no difference between the cumulative 
distributions of a v and <r v -s (using both the Kolmogorov- 
Smirnov and the Wilcoxon tests, e.g. Press et al. 1992), 
the median value of <j v /<t v -s being consistent with unity. 
Moreover, as for individual clusters, a v and a v s never 



differ by more than 10%. 

5 "ACTIVE" AND NON "ACTIVE" GALAXIES 

There is evidence that the spatial distribution and kine- 
matics of late-type galaxies (or blue galaxies or ELGs) are 
different from those of early-type galaxies (or red galax- 
ies or NELGs); these differences lead to higher estimates 
of internal velocity dispersion and mass for clusters (e.g., 
Moss & Dickens 1977; Biviano et al. 1997; Mohr et al. 
1996; Carlberg et al. 1997a; Koranyi & Geller 2000). Bi- 
viano et al. (1997) suggested that the dynamical state of 
the ELGs, which are often spirals of late types or irreg- 
ulars, reflects the phase of galaxy infall rather than the 
virialized condition in the relaxed cluster core. Carlberg 
et al. (1997a) suggested that, although differing in their 
distributions, both blue and red galaxies are in dynamical 
equilibrium with clusters. 

We check if the populations of "active" and non "active" 
galaxies, AGs and NAGs, really differ in their kinematics. 
Since we often classify AGs on the base of the presence of 
emission lines in their spectra (but see Couch et al. 1998 
for the use of colors, too) this classification roughly corre- 
sponds to the one between ELGs and NELGs. 




R/Rvir 

Fig. 3. — The (normalized) line-of-sight velocity dispersion, 
a v (R), as a function of the (normalized) projected distance 
from the cluster center. The points represent data combined 
from all clusters and binned in equispatial intervals. Open and 
filled circles are obtained using the "active" galaxies, AGs, and 
the galaxies without strong signs of activity, NAGs, respec- 
tively. The normalizing quantities are computed combining 
both the AGs and NAGs. We give the robust estimates of 
velocity dispersion and the respective bootstrap errors. 

We consider the 43 out of 51 clusters for which spectral 
information is available, each cluster containing Nnag(— 
N m ) NAGs and N AG (= N g - N m ) AGs, cf. Table 2. Fig- 
ure 3 shows that the <r v profile of the AGs is generally 
higher than the profile of the NAGs, where the profiles 
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are obtained combining together galaxies of all 43 clusters 
and normalizing to the values of a v and R V i r obtained 
for clusters before the rejection of the AGs (otherwise 
the difference would be also larger). A two-dimensional 
Kolmogorov-Smirnov test (Fasano & Franceschini 1987) 
applied to the normalized velocities and distances found 
a difference larger than > 99% between the two galaxy 
populations. 

Only for 19 of the 43 clusters there are enough galaxies 
(Nag > 5 & n d Nnag > 5) to compute and compare the 
respective AG and NAG u v . Figure 4 shows the compar- 
ison between cumulative distributions of a v as computed 
considering the AGs or NAGs. The AG a v -distribution 
shows a tail at high a v , which, however, is not significant 
according to a Kolmogorov-Smirnov test, and only slightly 
significant at the 93% c.l. according to the Wilcoxon test 
(e.g., Press et al. 1992). To test for different means and 
dispersions of the AG and NAG populations in each in- 
dividual cluster, we apply the standard means-test and 
F-test (e.g., Press et al. 1992). We find evidence of 
a difference more significant than > 95% only for A851 
(o'ag = 1761 km s _1 ^ onag = 1067 km s _1 at the 
- 98% c.l.) and for CL1602+4304 (< V >ag= 268653 
km s- 1 ^< V >nag= 270349 km s" 1 at the - 96% c.l.). 
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Fig. 4. — The cumulative distributions of line-of-sight velocity 
dispersion computed using the "active" galaxies, AGs, and the 
galaxies without strong signs of activity, NAGs, as indicated 
by dotted and solid lines, respectively. 

The evidence of difference between the two populations 
is in agreement with previous findings for nearby and dis- 
tant clusters (e.g., Biviano et al. 1997; Dressier et al. 
1999; Mohr et al. 1996). A really quantitative compari- 
son of the effect is complicated by the differences in the 
definition of "active" galaxies. By using ENACS clus- 
ters, Biviano et al. (1997) found that a v of ELGs is, 
on average, 20% larger than that of NELGs, and we find 
< 'AG I i<j n AG >= 1-12 ± 0.07. As for the combined veloc- 
ity dispersion profile (cf. Figure 3), our result in central 



cluster regions (R <^ 0.1 h 1 Mpc ) is similar to that of of 
Biviano et al (1997) giving higher central a v for AGs (and 
ELGs) than for NAGs (and NELGs), i.e. ~ 1.4 ± 0.2) 
vs. ~ 1.1 ± 0.05 for normalized a v . Our last point of 
the velocity dispersion profile is instead very high, but we 
suspect that it could be due to the loss of efficiency in re- 
jecting interlopers in very poorly sampled external regions 
of distant clusters. 

6 COMPARISON WITH RESULTS 
FROM X-RAY AND LENSING DATA 

We collect X-ray luminosities, in general bolometric 
ones, Ltoi.x, and temperature, Tx, for 38 and 22 clus- 
ters, respectively. For A223, A521, CL0054-27, CL0412- 
65, CL1604+4304, and the four clusters by Bower et al. 
(1997) we obtain the bolometric luminosities by multi- 
plying the original band luminosities by a temperature- 
dependent bolometric correction factor. This factor is 
computed under the assumptions of pure bremsstrahlung 
intraclustcr medium emission and a power-law approxi- 
mation for the Gaunt factor. For the correction we use 
the temperatures estimated from a v in the hypothesis of 
density energy equipartition between hot gas and galax- 
ies, i.e. (3 ' S p ec — al / (kT / fim p ) = 1, where fj, = 0.58 is the 
mean molecular weight and m p the proton mass. For the 
four clusters by Bower et al. (1997), which have very few 
galaxies with measured redshift, we use the Lx — cr v re- 
lation given by the authors for nearby clusters. For these 
four clusters we expect a v ~ 600 km s _1 , i.e. T ~ 2 keV. 

In Table 4 we list the collected values for L^lx, and 
Tx with the corresponding reference sources (Cols. 2-5), 
and the value of p spec (Col. 6). The errors on I3 spec take 
into account errors on both a v and Tx- 

Figure 5 shows the Lx,boi-&v relation compared to that 
found by Borgani et al. (1999) for nearby clusters. Ex- 
cluding the leftmost point (J2175.15TR), the resulting bi- 
secting linear regression is 



-i2.6±i-.2, 



logiLbd.x/W^ergs- 1 ) = 4A±\ilog(a v /km s" 1 )- 

(6) 

where errors come from the difference with respect to the 
direct and the inverse linear regression (Isobc et al. 1990, 
OLS methods). Our Lf, i,x~&v relation is consistent with 
that of nearby clusters (e.g., White et al. 1997; Borgani 
et al. 1999; Wu, Xue, & Fang 1999). As for the point 
excluded, note that our analysis of J2175.15TR is based 
only on 19 galaxies, and the estimate of a v is recovered 
from only eight member galaxies (with an error larger than 
100%). 

All clusters for which a v should be better inter- 
preted as an upper limit to the true estimate (AS506, 
CL0054-27, 3C295) lie on the right-upper corner of the 
plot. Excluding these points we fit a consistent relation, 
i.e. logiL^x/lO^ergs- 1 ) = 4.7t{ilog(a v /km s' 1 ) - 
13.5±|"i. 

Figure 6 shows the a v -Tx relation compared to that of 
nearby clusters, as reported by G98. As for distant 
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TABLE 4 

Comparison with X-ray Properties 
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a Bolometric luminosity here computed from the original band luminosity. 

a Errors on Tx are at the 90% c.l.; they are rescaled by a factor 1.6 to compute 68% c.l. errors on /3spec ■ 

(1) White ct al. 1997; (2) Wu ct al. 1999 (see this compilation for original data sources); (3) Lea & Henry 1988; (4) Ulmer et al. 1985; (5) Soucail 
et al. 2000; (6) Smail et al. 1997; (7) Kaiser ct al. 1998; (8) Castandcr ct al. 1994; (9) Tucker et al. 1998; (10) Bower ct al. 1997; (11) Pierre et 
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clusters, the data have a too small dynamical range to 
attempt a linear fit: the visual inspection of Figure 6 sug- 
gests no difference with nearby clusters in agreement with 
the result of by Mushotzky & Scharf (1997) and Wu ct 
al. (1999). We obtain f3 spec = 0.88+^17 (median value 
with errors at the 90% c.l.), in agreement with the value 
of (i spec = 0.88 ±0.04 for nearby clusters (cf. G98). More- 
over, we find no correlation between fi spec and redshift (cf. 
also Wu et al. 1999). 

Under the assumption that the hot diffuse gas is in hy- 



drostatic and isothermal equilibrium with the underlying 
gravitational potentials of clusters, one can obtain the X- 
ray cluster masses provided that the gas temperature and 
radial profile of gas distribution are known (e.g. Wu et 
al. 1998). The availability of Tx allow us to compute 
the mass within R V i r for 22 clusters according to Mx = 

where we adopt the gas distribution given by the /3-modcl 
with typical parameters (slope Pfa,gas = 2/3 and core ra- 
dius R x = 0.125 h,- 1 Mpc , e.g. Jones & Forman 1992). 
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We find mass values consistent with our optical virial esti- 
mates, i.e. M/M x = 1.02 (0.86-1.32) for the median value 
and the range at the 90% c.L. 



100 (= r 




1000 

ct v / (km s" 1 ) 

Fig. 5. — L x ,boi-o~v relation for distant (open circles) and 
nearby clusters (filled circles) . For the nearby clusters we show 
results as reported by Borgani et al. (1999), all having a v esti- 
mated at least with 30 galaxy redshifts (Girardi et al. 1998b) 
and also belonging to the X-ray Brightest Abell-like Cluster 
survey (Ebeling et al. 1996). The error bands at the 68% c.L 
are shown: errors on Lx,boi are not available for a few distant 
clusters. The three solid lines are direct, inverse, and bisecting 
linear regression for the distant clusters (obtained rejecting the 
point on the left). The dashed line is the bisecting linear re- 
gression for the nearby clusters as computed by Borgani et al. 
(1999). 

As for gravitational lensing masses, we resort to esti- 
mates found in the literature. We collect projected esti- 
mates from weak gravitational lensing analysis, Ml, for 
18 clusters. In order to compare our optical virial masses 
to Ml, we project and rescale our masses M within the 
corresponding radius using the fitted galaxy spatial distri- 
bution. In Table 5 we list the reference sources (Col. 2) 
from where we take Rl and Ml (Cols. 3 and 4); the cor- 
responding optical virial projected mass, M opt ,L (Col. 5); 
and the respective ratio (Col. 6). We obtain M opt .L/ML = 
1.30(0.63-2.13) (median value and range at the 90% c.L). 
Moreover, we do not find any correlation between M/Mx 
or M optt LlML and redshift. 

Our finding are in agreement with other recent stud- 
ies which find, on average, no evidence of discrepancy be- 
tween different mass estimates as computed within large 
radii, thus suggesting that distant clusters are nor far from 
global dynamical equilibrium (e.g., Allen 1998; Wu et al. 
1998; Lewis et al. 1999). Note that we avoid to consider 



mass determination in very central cluster regions since 
our analysis of cluster members give poor constrains on 
mass distribution on these scales. Indeed, the assumption 
of dynamical equilibrium seem to break down in the very 
central regions as suggested by comparisons with strong 
lensing mass estimates (e.g., Allen 1998; Lewis et al. 1999; 
Wu et al. 1998). 




1 10 

kT x / keV 

Fig. 6. — a v -Tx relation for distant (open circles) and nearby 
clusters (filled circles) . For the nearby clusters we show results 
as reported by Girardi et al. (1998b), all having o v estimated 
at least with 30 galaxy redshifts, and Tx taken from David et 
al. (1993) and from White et al. (1997). The error bands at 
the 68% c.L are shown: when authors give only 90% c.L errors 
on Tx, we apply a reduction by a factor of 1.6. The solid line 
is the bisecting linear regression for the nearby clusters as com- 
puted by Girardi et al. (1998b). The dashed line represents the 
model with the equipartition of energy per unit mass between 
gas and galaxy components (/3 spec = 1). 



7 SUMMARY AND CONCLUSIONS 

In order to properly analyze the possible dynamical evo- 
lution of galaxy clusters we apply the same procedures al- 
ready applied on a sample of nearby clusters (170 clusters 
at z < 0.15 from ENACS and other literature, Girardi et 
al. 1998b, cf. also Fadda et al. 1996) to a corresponding 
sample of 51 distant clusters (0.15 <; z <; 0.9, < z >^ 0.3). 
Each cluster has at least 10 galaxies with available redshift 
in the literature. Three cluster fields show two overlapping 
peaks in their velocity distribution and large uncertainties 
in their dynamics. Out of other cluster fields, 45 fields 
show only one peak in the velocity distribution and three 
fields show two separable peaks, for a total of 51 well- 
defined cluster systems. These 51 systems are those 
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TABLE 5 

Comparison with Masses from Weak Gravitational Lensing 
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(1) Seitz et al. 1996; (2) Squires et al. 1996a; (3) Squires et al. 1996b; (4) Natarajan ct al. 1998; (5) Smail ct al. 1997; (6) Kaiser et al. 1998; 
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used in the comparison with nearby clusters (i.e., 160 well- 
defined systems, cf. § 3 of Girardi et al. 1998b). 

We select member galaxies, analyze the velocity disper- 
sion profiles, and evaluate in a homogeneous way cluster 
velocity dispersions and virial masses. 

As a main general result, we do not find any significant 
evidence for dynamical evolution of galaxy clusters. More 
in detail, our results can be summarized as follows. 

• The galaxy spatial distribution is similar to that of 
nearby clusters, i.e. the fit to a King-like profile 
gives a two-dimensional slope of a — 0.7 and a very 
small core radius of R c = 0.05 h^ 1 Mpc. Note that 
we do not want to really asses the existence of a 
core, or to state that the King-modified profile is 
better than other forms for galaxy density profiles; 
the King-modified profile is used for a consistent 
comparison with nearby clusters. We refer to Gi- 
rardi ct al. (1998b, § 8) for other relevant analyses 
and discussions. 

• For those clusters with good enough data, the inte- 
grated velocity dispersion profiles show a behavior 
similar to those of nearby clusters: they are strongly 
increasing or decreasing in the central cluster re- 
gions, but always flattening out in the external re- 
gions, thus suggesting that large-scale dynamics is 
not affected by velocity anisotropies. 

• The average velocity dispersion profile can be ex- 
plained by a model with isotropic orbits, which well 
describe also nearby clusters. Possible evidences for 
more radial orbits are not statistically significant. 



• There is no evidence of evolution in both Lbo^x-Cv 
and a v -Tx relations, thus in agreement with previ- 
ous results (Mushotzky & Scharf 1997; Borgani et 
al. 1999). 



Moreover, within the large scatter of present data, we 
find, on average, no significant evidence of discrepancies 
between our virial mass estimates and those from X-ray 
and gravitational lensing data, thus suggesting that dis- 
tant clusters are not far from global dynamical equilibrium 
(cf. also Allen 1998; Lewis et al. 1999; Wu et al. 1998). 

We conclude that the typical rcdshift of cluster for- 
mation is higher than that of our sample in agreement 
with previous suggestions (e.g., Schindler 1999; Mushotzky 
2000). In particular, we agree with preliminary results by 
Adami et al. (1999), who applied the same techniques 
used for the nearby ENACS clusters on 15 distant clus- 
ters, (< z >^ 0.4) from the Palomar Distant Cluster Sur- 
vey (Postman et al. 1996). 

Although some clusters at very high redshift, e.g. z > 
0.8, are already known (e.g., Gioia et al. 1999; Rosati et 
al. 1999), the construction of a large cluster sample useful 
for studying internal dynamics will require a strong ob- 
servational effort. Note that, already in the construction 
of the cluster sample analyzed here, we relax the require- 
ments applied to the sample of nearby clusters by Girardi 
et al. (1998b), i.e. the distant clusters are more poorly 
sampled. Throughout the presentation of our analyses we 
stress how both the poor number of galaxies and the small 
spatial extension of some clusters can affect the robustness 
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of their resulting properties. In particular, Monte Carlo 
simulations, which take into account the whole member- 
ship procedure, show that the estimate of velocity disper- 
sion is, on average, well recovered also in the case of very 
poor sampling (only 10 galaxies in the cluster field giving 
5-6 members), but that the global error associated to the 
individual clusters should be a factor ~ 2.5 larger than 
the pure statistical error. Also the small spatial extension 
could lead to large over/underestimates of velocity disper- 
sion of individual clusters: we evaluate that variations of 
25% are quite common when clusters are sampled only out 
to half of the virialization region. 
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with the National Aeronautics and Space Administration. 
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Ministry of University, Scientific Technological Research 
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Appendix A 
RESULTS FOR 
MULTIPEAKED CLUSTERS 

Here we shortly present the results of our analysis for the 
three clusters with uncertain dynamics, i.e. with two peaks 
in the velocity distribution which are not clearly separable. 
We consider both the system composed by the two peaks 
together and each peak individually. In Figure 7 we plot 
the velocity-space galaxy density for each cluster, as pro- 



vided by the adaptive-kernel reconstruction method, and 
the integrated velocity dispersion profile VDP for each pos- 
sible system. Table 6 summarizes the results of the analy- 
sis of the internal dynamics. Note the strong variation in 
a v and mass when considering the two peaks together or 
not. Some comments on individual clusters follow. 
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Fig. 7. — For each of the three clusters we give the relative 
velocity-space galaxy density, as provided by the adaptive ker- 
nel reconstruction method (left panels) and integrated line-of- 
sight velocity dispersion profiles a v (< R) for each of the consid- 
ered systems (right panels). For the velocity dispersion profiles 
we plot bootstrap errors only in the case of the system with 
joined peaks. The horizontal lines in the right panels represent 
X-ray temperature with the respective errors taken from Wu 
et al. (1999) and transformed in a v imposing f3 spec = 1. 



TABLE 6 

Clusters with Uncertain Dynamics 
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A1689 Teague, Carter, & Gray (1990) computed a value 
of <r v = 1989 km s _1 . As for the analysis of the cluster 



members, the two peaks in the velocity distribution were 
already pointed out by Girardi et al. (1996; 1997b) us- 
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ing the same adaptive kernel method. By using a multi- 
scale analysis which couples kinematical estimators with 
the wavelet transform, Girardi et al. (1997b) found the 
presence of two dynamically relevant structures, but with 
a smaller a v and mass with respect to the two systems, "a" 
and "b", analyzed here. Moreover, A1689 is well know 
for a strong discrepancy between mass from X-ray and 
strong gravitational lensing analyses (e.g. Miralda-Escude 
& Babul 1995) which could be due to its complex struc- 
ture. Also the very recent weak lensing analysis of Taylor 
et al. (1998) suggests the model of a double cluster aligned 
along the line of sight in order to explain discrepancies 
between optical and X-ray results. These results and the 



fact that A1689 appears well aligned along the line of sight 
with other structures (three foreground groups, Teague et 
al. 1990) suggests the presence of a large structure fila- 
ment well aligned along the line of sight. 
A2744 Couch & Sharpies (1987) computed a value of 
<j v = 1947^201 km s _1 . A strong suggestion for the dy- 
namical activity comes from the recent study by Allen 
(1998): among the 13 clusters analyzed, A2744 shows the 
strongest discrepancy between mass from X-ray and grav- 
itational lensing analyses. 

A3854 Colless & Hewett (1987) listed a value of a v = 
1180;??o km a" 1 . 
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